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We report on the largest Swift AGN monitoring program, concerning UV/optical variability in 
Seyferts. From 554 observations, over a 750d period, of the Seyfert galaxy NGC 5548, we see 
(McHardy et al. 2014) a good overall correlation between the X-ray and UV/optical bands, par¬ 
ticularly on short timescales (tens of days). The UVOT bands are found to lag behind X-rays 
with a lag scaling as wavelength to the power 1.23 +/- 0.31, in excellent agreement with that ex¬ 
pected (1.33) if UV/optical variability arises from reprocessing of X-rays by the accretion disc. 
However, the observed lags are ~ 3 times longer than expected from a standard Shakura-Sunyaev 
disc, raising real concerns about the detailed validity of this model. The results can be explained 
with a slightly larger mass and accretion rate, and a hotter disc, or if the disc is clumpy, thereby 
enhancing the emission from the outer regions. 
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1. Introduction 

The relationship between UV/optical variability and X-ray variability in AGN, and the origin 
of the UV/optical variability, is still not clear. Strong X-ray/UV or X-ray /optical correlations, with 
short lags of less than a day, have been observed on short timescales (weeks - months) by a number 
of studies [e.g. 24, 2, 4, 15, 7, 8, 22], However, these studies show poorer correlations on longer 
timescales (months - years), usually due to long-term UV/optical trends not present in the X-ray 
variability, suggesting that UV/optical variability is dominated by different processes on different 
timescales. Reprocessing of X-rays into the UV/optical emission by the nearby accretion disc could 
cause an X-ray/UV correlation on short timescales, leading to UV/optical variations lagging the X- 
rays by the short light travel time between the X-ray source and the disc. Alternatively, if X-ray 
variability is produced by variability in the UV seed photon flux due to accretion variations in the 
disc at very small radii, a correlation could arise in which X-ray variability lags the UV-optical, by 
that same light travel time. Determining the precise lag between variations in the X-ray and UV- 
optical emission is therefore extremely useful in determining the origin of UV-optical variability. 



Figure 1: (Bottom panel) Long term Swift 0.5-10 keV X-ray count rate. (Top Panel) UVW2 flux. 

Short (~ 1 d) lags of the X-rays by the optical are seen in all previous studies, which mostly 
consist of a combination of RXTE X-ray monitoring and ground-based optical monitoring. The lag 
has, however, never been measured well enough to completely rule out an optical lead unambigu¬ 
ously. Whereas ground-based monitoring is affected by weather, the Swift observatory can provide 
uninterrupted simultaneous X-ray and UV/optical monitoring, allowing more precise measurement 
of wavelength-dependent lags. 

Using Swift observations, the B-band was shown to lag the X-rays by < 45mins in NGC 4395 
[7] and interband lags were measured in NGC 2617 [22] which agree well with predictions from a 
reprocessing model. 

Here, we discuss the largest Swift AGN monitoring campaign to date, of the Seyfert 1 galaxy 
NGC 5548, consisting of 554 observations over a 750d period [17]. The observations were not 
scheduled to follow particular events, but are instead typical of the long-term behaviour of NGC 
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Figure 2: Multiband UVOT light curves in mJy. 

5548. A strong X-ray/V-band correlation was already known to exist in NGC 5548 [24], but until 
this work the V-band lag had not yet been precisely defined. 


2. Observations & Data Reduction 


The SWIFT X-ray data were taken by the X-ray Telescope (XRT, [5]) and the UV/optical 
observations taken by the UV and Optical Telescope (UVOT, [20]). All XRT observations were 
carried out in photon-counting (PC) mode and all UVOT observations were earned out in image 
mode. The data were analysed using our automatic pipeline, based upon the standard Swift analysis 
tasks, as described in e.g [7, 10]. The X-ray data are corrected for any effects caused by vignetting 
and aperture losses and bad pixels. Drops in UVOT data points were investigated individually, and 
those which were due to bad tracking or bad pixels were removed. 

Observations took place between MJD-50000 of 5960-6709, typically every 2 days, with some 
periods of less or more frequent sampling (Id or 4d), and were mostly of 1 ks though sometimes of 
2 ks. Observations were usually split into 2 or more individual visits, improving time sampling. A 
total of 554 visits were made, giving, after rejection of bad data, 465 usable X-ray measurements. 

As found previously between the X-ray and V bands [24], we find fairly close correspondence 
between the X-ray and UVW2 flux lightcurves. On short timescales (~ lOd), the correspondence 
is strong, but on longer timescales the amplitudes of variability are not always equal. Observations 
were made in additional UVOT filters from Day 6380; the resultant light curves are shown in Fig. |j. 
A very close correspondence is apparent between all UVOT bands. 


3. X-ray / UV-Optical Correlations 

3.1 The X-ray/UVW2 lag 

In Fig. |3] we show the discrete [DCF, 12] and interpolation [ICCF, 14] CCFs of the UVW2 
band, which is the best sampled of the UVOT bands, for the complete X-ray and UVW2 datasets. 
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Figure 3: Left: Discrete cross correlation function between the X-ray and UVW2 lightcurves shown in 
Fig. 0 - The 95% (dashed red) and 99.99% (solid thin blue) confidence levels are also shown. Right: In¬ 
terpolation cross correlation function between the X-ray and UVW2 lightcurves shown in Fig. [I|. The 95% 
(dashed red) and 99% (dashed blue) confidence levels are also shown. 



Lag of X-rays by UVW2 (days) Lag of UVW2 by V (days) 


Figure 4: Lag distributions from Javelin: Left Panel UVW2 following X-rays using data, following mean 
subtraction, from the intensive period from 6383-6547. Right Panel V following UVW2 using data shown 
in Fig. | 


Only the UVW2 band, being more slowly varying, is interpolated. The weighted mean ICCF of 
the 3 sections seen in Fig. [l] is taken, rather than interpolating across the gaps. The N% confidence 
curves are also shown, and are defined such that (100-N)% of correlations performed between the 
UVW2 data and randomly simulated X-ray lightcurves with the same variability properties as the 
X-ray data would exceed these levels [e.g. 23], (e.g. see [4] for more details). Both the DCF and 
ICCF show a broad, but highly significant correlation, with a peak near a lag of zero, with the DCF 
favouring a lag of the X-rays by the UVW2 of about a day. 

As correlation functions can be distorted by a long term variations in the mean level not present 
in both lightcurves, it is recommended practice to subtract a running mean [26]. A mean based on 
a running boxcar of width 20d is therefore subtracted from both UVW2 and X-ray light curves. 

In order to refine the lag measurement, we calculated the distribution of lags from 10000 
simulations for the period of intensive observation (Day 6383 - 6547), using the JAVELIN cross¬ 
correlation program [28], used previously in e.g. [22, 19], allowing a lag range of ±10d. JAVELIN 
assumes variability based on a damped random walk in both bands [27]; although the X-ray and 
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Figure 5: Lag of the X-ray and other UVOT bands relative to UVW2. Lag oc wavelength^ where (i = 
1.23 ±0.31. The lower, dashed, red line is the prediction for a standard disc as described in the text. 

UVW2 variability properties are slightly different, lag measurements are not significantly affected 
[22]. The resultant lag distribution, shown in Fig. ^ (left panel), has a median UVW2 lag of 
+0.70^02^. When applied to the non-mean subtracted light curves, JAVELIN does not converge 
to the same single distribution, presumably due to the uncorrelated long timescale variations. 

3.2 UVOT interband lags 

As the UVOT bands all show the same long term trends, it was not necessary to mean subtract 
the lightcurves. The lag distributions between UVW2 and the other UVOT bands were therefore 
calculated with JAVELIN using data shown in Fig. |2[ an example of which is shown in Fig. [| 
right panel. As some distributions, are slightly asymmetric, use of the mode rather than the median 
would slightly reduce the lags in some cases, but the differences are small. 

3.3 X-ray reprocessing in a ‘Standard Disc’ 

On average, the UVOT bands lag the X-ray band with lags which increase with wavelength 
(Fig.|). These lags can be compared to predictions of reprocessing from a simple accretion disc, for 
which the lag should vary as the 4/3 power of wavelength [e.g. 6, 15]. The lags were therefore fitted 
with a simple model of the form lag = A + (B x A )P. The best fit, assuming Gaussian distributed 
errors, gives parameters of A = —0.70 ± 0.21, /3 = 1.23 ± 0.31 and B = 3.2 ±0.6 x 10~ 3 , with a 
X 2 of 2.1 for 3 d.o.f. The slope, /3, is similar to that derived by [22] for NGC 2617 (1.18 ±0.33), 
but does not require an unphysical offset to the X-ray point, as they required. 

This result demonstrates that the short term UV/optical variability in NGC5548 can be very 
well explained by reprocessing of X-ray emission from an accretion disc. 

4. Accretion disc Modelling 

In order to determine whether the observed lags can be explained by a disc consistent with 
other observed properties of NGC 5548, we model the disc using the methods described in [15] 
(c.f. [8]), including the effects of both X-ray and gravitational heating. The key parameters in the 
model are the X-ray heating, obtained by extrapolating the observed 2-10 keV luminosity (Lx 2 - 10 ) 
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to ~0.01-500 keV, the disc albedo, the height, H , of the X-ray source above the disc (assuming a 
lamppost geometry), and the inner disc radius, R m . 

We assume a black hole mass, Mbh, of 6.7 x 10 7 M 0 [3], an accretion rate, the, of ~ 0.03 — 
0.04 of Eddington [25], a high X-ray heating luminosity 6 x Lx 2-10 (implying a low albedo of 
20%), H = 6 7C,(consistent with X-ray source sizes measured by other methods [9, 13]) and R,,, = 6 
R g , the ISCO for a Schwarzschild black hole. This produces the dashed line shown in Fig. |j. The 
lags predicted from this standard disc model represent the time for half of the reprocessed light has 
been received, following impulse X-ray illumination. As the response is asymmetric, however, the 
peak response may be even faster, giving even shorter predicted model lags. 

In order to increase the predicted lags from this homogeneous disc model such that they agree 
with observation, the geometry would have to be changed (e.g. H = 20 R g , R nl = 20 R g ), and a 
larger Mbh = 10 8 M 0 ) and hotter (m/j = 0.06) disc are required. In this case, disc flaring would 
have a smaller effect than an increased disc temperature. 

5. Discussion 

5.1 Long X-ray/UV-Optical Lags 

The form of the observed wavelength-dependent lags, i.e. lag ~ A 1 23 (Fig. |||), strongly support 
reprocessing of X-rays as the main cause of short-timescale UV/optical variability in AGN, how¬ 
ever our observations also imply that this reprocessed emission originates further from the black 
hole than predicted by a ‘standard’ accretion disc model. The predictions and observations can 
only be reconciled by pushing parameter limits. It is important to note, however, that a larger than 
expected disc is also required [18] to explain microlensing observations, lending further support to 
the idea that the standard disc model is often inadequate, as proposed by e.g. [1]. The suggestion 
[18] of a low radiative efficiency is consistent with our requirement for R, n > 20. Alternatively, an 
inhomogeneous discs, for which the outer portions contribute more flux than for a uniform disc, 
would cause the disc to appear larger [11]. 

5.2 Long-Timescale UVOT variability 

Although the correspondence of X-ray/UVW2 is very good on short timescales, and generally 
good on long timescales, from Day ~ 6470 to 6547 there is a rise in the UV/optical which is sig¬ 
nificantly more pronounced than in X-rays. The UV/optical rise could be explained by an inwardly 
propagating rise in accretion rate, which eventually hits the X-ray emission region, however, the 
viscous timescales from 100 R g are too long (~ 10 years, c.f.[4]). The data can therefore only be 
explained if the ratio of disc scale height to radius is larger than the value of 0.1 which is normally 
assumed, or the rise propagates through a corona over the disc rather than the disc itself, or the 
perturbation starts at a small radius (~ 20 R g ). 

5.3 Seed Photon Variations as a Driver of Variability 

Although the result that the X-rays always lead the UV/optical can be explained by assuming 
that all X-ray variability is generated within the corona, it is still interesting to ask why variations 
in seed photon flux, which would cause the UV-optical emission to lead, appear to have little affect 
on the measured lags. We suggest that the answer is a combination of relative solid angles and 
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conservation of photons during the X-ray scattering process. A given single seed photon, Compton 
up-scattered into an X-ray photon with energies 10-100 x greater than the seed photon, can, through 
reprocessing in the disc, produce many UV-optical photons by black body emission. As the disc fills 
half the solid angle seen by the X-ray source, a large fraction of the X-ray photons will hit the disc. 
The total number of UV-optical photons produced in this way may therefore exceed the initial seed 
photon fluctuation. The reprocessing time, which is usually ignored in reprocessing calculations 
but which must be of finite length, could also add a further delay to the optical lightcurves, aiding 
agreement with observation. 
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